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We report results of fully non-perturbative calculations, based on Auxiliary Field Quan- 
tum Monte Carlo (AFQMC) approach, for the dilute neutron matter at the density 
p = 0.003 fm~ 3 . Fundamental quantities which characterize the superfluid state: the 
single particle energy gap A(T), and the critical temperature T c have been determined. 
The large value of A(0)/T c ?s 3.2 indicates that the system is not a BCS-type superfluid 
at low temperatures. 

1. Introduction 

Although a homogeneous neutron matter is one of the simplest nuclear systems, its 
importance cannot be overemphasized since it constitutes the main component of 
neutron stars. The density of the neutron matter forming neutron stars ranges from 
subnuclear densities (in the inner crust) up to extremely high values of p ~ 10 on 
(expected in the center of the star), where pa = 0.16 fm -3 is the saturation density-M 
In the regime of sufficiently low densities one can perform very precise cal- 
culations since the neutron-neutron interaction is completely dominated by the 
scattering in So channel, which is determined by two parameters only: the scat- 
tering length a and the effective range r ff. Indeed, at the densities p < 0.03 po 
(PF i5 0.6 fm -1 ) the influence of other channels as well as of three-body forces is 
marginal and can be neglected^. The values of the scattering length and the effec- 
tive range are well known from the low energy scattering experiments and read: 
a = -18.8 ± 0.3 fm, r cff = 2.75 ± 0.11 fm 3 . At the densities where |p F a| > 1 the 
dilute neutron matter is an example of a strongly correlated system. Moreover, 
the influence of the effective range cannot be ignored except for very low densities 
where pF^eff << 1 ■ It implies that only non-perturbative approaches are able to 
gain a reliable insight into physics of this system. The large class of such methods, 
which are known under the general name of Q uantum Monte Carlo (QMC), has 
been applied for the dilute neutron matter^EEl i n spite of a considerable theoret- 
ical effort some of its static properties are still vaguely known. The open questions 
mostly concern the superfluid state, which is generated at sufficiently low temper- 
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atures. For example the predicted value of the zero temperature energy gap differs 
significantly for various approaches, revealing a strong dependence on the details 
of the applied method (see^ and references therein) . The same concerns the critical 
temperature of the superfluid-normal phase transition, where the well known BCS 
formula A(0)/T c w 1.76 is believed to be valid. 

Here we present the selected results concerning the superfluid state, obtained 
within the AFQMC approach at non-zero temperatures^^ for the dilute neutron 
matter of density p — 0.02 po (pp = 0.45 fm _1 ). This study can be regarded as a 
continuation of the investigations presented in the papers 021 where the details of 
the AFQMC algorithm and the procedure of adjusting the interaction, have been 
discussed. In the present study we focused on two basic quantities which characterize 
a superfluid state: the critical temperature and the pairing gap. 

2. Critical temperature 

In the case of superfluidity the convenient order parameter is specified by 
the long-distance behavior of the two-body density matrix P2{^i, ri, rj) = 
(V'l(ri)^l (13)^4.(^2)^ ( r i))' wnere the field operators ifi\{r) and ip\(r) obey the 
fcrmionic anticommutation relations. Such a long-distance behavior is directly re- 
lated to the condensate fraction a, which measure the fraction of Cooper pairs 
forming a condensate^. This quantity can be efficiently computed within AFQMC, 
provided the information concerning the one-body density matri x O ] (ri , ra) = 
(i/'j(r2)^4,(ri)} is also included. The condensate fraction is given byP^ 



where N stands for the number of particles and the correlation functions g-yp ( r ) 
arc defined as: 



The figure [T] presents the condensate fraction versus the temperature for three 
different lattice sizes N s . It is clearly visible that the results for N s = 6 significantly 
deviate from those obtained for larger lattices. It is due to the fact that for small 
lattices the limiting value of h(r) is rather poorly determined. Indeed, "the infinity" 
in this case is equal to the half of the box size and is only a few times larger than the 
interaction range. Consequently, it will generate a large systematic error when the 
lattice is too small. The presented results indicate that this is exactly the case for 
N s = 6 and therefore these data were not included in the process of determination 
of the critical temperature. 

In the thermodynamic limit the condensate fraction has to vanish at the critical 
temperature T c . However, in the case of calculations in the box the finite size effects 



a = lim h(r), h{r) = — g 2 (r) - gf(r), 
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Fig. 1. The condensate fraction a as a function of temperature (in units of Fermi energy ep) 
for three lattice sizes N s . The inset presents the rescaled condensate fraction obtained for lattices 
N s = 8 (squares) and N B = 10 (circles). The intersection of the fitted lines determines the critical 
temperature of the supcrnuid-normal phase transition, T c S3 0.09 £f. 

come into play and smooth out all singularities typical for the phase transition. Still, 
one can estimate the critical temperature using the method based on the finite size 
scaling. Similar techniq ue h as been used to determine the critical temperature for 
cold atomic gases for details). The volume-dependent estimation of the 

critical temperature is obtained by finding the intersection of the rescaled 

condensate fraction L 1+n a for two different lattice sizes Nij, where r] ~ 0.038 is 
the universal critical exponent and L 3 is the volume of the box. As Nij — > oo the 
thermodynamic limit is recovered and the series converges to the critical tem- 
perature. To extract T c we have used two largest of the available lattices, namely 
N id = 8, 10. The filling factor in both cases reads: v = N/2N 3 w 5%. Accord- 
ing to our experience, the above methodology is sufficient to estimate the critical 
temperature with an uncertainty smaller than 20% (in fact this procedure applied 
to the unitary gas provides the value of T c with the relative error smaller than 
10%). Finally, the estimated value of the critical temperature reads T c « 0.09 £f 
(0.39 MeV), where Ef is the Fermi energy (see the inset of Fig.[T]). 

3. Pairing gap 

In order to determine the gap in the single particle spectrum, we have computed the 
spectral weight function A(p,uj). This quantity has been extracted from the imag- 
inary time propagator Q{p,r) through the analytic continuation^^!. The procedure 
is equivalent to solving the integral equation: 

2tt J_ oa 1 + exp(-cj/3) 



i 
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(w-p)/e F {u-n)/e F 



Fig. 2. The spectral weight function yl(p,a;) for the dilute neutron matter at selected tempera- 
tures: A) T ~ 0.06e P < T c , B) T ~ 0.08e F § T c , C) T ~ 0.10e F § T c , D) T ~ 0.12e F > T c . The 
nonzero energy gap is present up to the critical temperature T c 0.09 ep. 



where j3 denotes the inverse of the temperature and G(p,t) is determined from 
the Monte Carlo calculations for the discrete set of values To = 0,n, . . . , r„ T = (3. 
Numerically, however, the above integral equation represents an ill-posed problem 
i.e. there is an infinite class of solutions for A(p,uj) which satisfy Eq. Q within 
uncertainties generated by Monte Carlo method. Therefore we have used two meth- 
ods which were in particular designed to d eal with such tasks: the truncated SVD 
method and the maximum entropy 

met hocp™a 

The spectral functions for selected temperatures (below and above T c ) are pre- 
sented in Fig. [21 Our results admit the gapped solutions up to the critical tempera- 
ture and above T c the gap vanishes. One has to remember, however, that there is a 
finite resolution of both methods concerning the value of the gap. Namely, the gap 
cannot be resolved if its value decreases below A min « 0.2 ep. 

The value of the energy gap A extracted from the spectral function at the tem- 
perature T ~ 0.06 £f is expected to provide a good approximation of its value at 
zero temperature. It reads: A(0)/ef = 0.29^q 04- ^ n Fig.[3]our results are compared 
to those obtained by other authors within ab initio calculations for zero tempera- 
ture. We have found our results in agreement with the recent Green Function Monte 
Carlo calculations of Gezerlis and Carlson^. The calculations including other scat- 
tering channels together with the three body forces were performed by Gandolfi et 
aL,^and predict larger paring gap. However, it is difficult to establish the source 
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Fig. 3. Supcrfluid paring gap in various Monte Carlo calculations together with results of pure BCS 
approach (dashed line) - see text for details. The dotted line corresponds to A(pp) = 0.29pp/2m. 

of this discrepancy since their method (Auxiliary Field Diffusion Monte Carlo) has 
been constrained in order to avoid the fermionic sign problem and therefore can be 
regarded as a variational approach. 

An interesting result can be noticed by considering the ratio of the energy gap 
at zero temperature to the critical temperature. Namely, the ratio A/T c w 3.2 
exceeds almost twice the well-known BCS value 1.76. A similar situation has been 
encountered for the unitary gas, where the existence of the exotic "pseudogap" 
phase above T c was recently reported^. It indicates that the dilute neutron matter 
at this density is not a BCS-type superfluid. 

4. Conclusions 

We have performed the fully non-perturbative calculations for the dilute neutron 
matter of density p = 0.02po at finite temperatures. We focused on the basic quan- 
tities which characterize the superfluid state and extracted the value of the critical 
temperature and the pairing gap. Our results are free from uncontrolled approxima- 
tions and are essentially exact with only uncertainties related to statistical errors 
and finite size effects. The large value of A/T c w 3.2 suggests that the dilute neu- 
tron gas cannot be described by the BCS theory and thus is not a BCS superfluid 
at low temperatures. 
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